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Abstract 

We consider the dynamics of a quantum particle in a one-dimen- 
sional periodic potential (lattice) under the action of a static and 
time-periodic field. The analysis is based on a nearest-neighbor tight- 
binding model which allows a convenient closed form description of 
the transport properties in terms of generalized Bessel functions. The 
case of bichromatic driving is analyzed in detail and the intricate trans- 
port and localization phenomena depending on the communicability of 
the two excitation frequencies and the Bloch frequency are discussed. 
The case of polychromatic driving is also discussed, in particular for 
flipped static fields, i.e. rectangular pulses, which can support an al- 
most dispersionless transport with a velocity independent of the field 
amplitude. 



1 Introduction 



Quantum transport properties in periodic structures (lattices) are highly non- 
intuitive, in particular in view of their localization properties (see, e.g., [1]). 
Classically, a particle in a periodic potential is localized in a single well at low 
energies, whereas a quantum system shows the well known Bloch band/gap 
structure with transporting bands. If a constant field Fq is added, the quan- 
tum states in a band with width G become localized to an interval G/Fq, 

*Present address: Institut fiir Physik, Universitat Kassel, D-34109 Kassel, Germany 



1 



contrary to our intuition. In this region, the wave functions show the cele- 
brated Bloch oscillations [2] with the Bloch frequency lob — F^d/h, where d 
is the period of the potential, and there is no directed transport. This local- 
ization is, of course, approximate because of the decay to infinity by Zener 
tunneling [3] , an effect of longer time scales at least for weak fields. An addi- 
tional time-periodic driving F{t) re-introduces the quantum transport which 
can be suppressed again for special choice of the parameters. In the following 
we will neglect the decay and confine ourselves to the minimal model system, 
the single-band tight-binding hamiltonian with nearest-neighbor coupling 

„ +00 +00 

i=—oo £=—00 

where \i) are the Wannier states which are exponentially localized in the 
i-th potential well and G is the width of the Bloch band for F = 0. Let us 
recall that for a static field F = Fq the extension of the Bloch oscillation is 
approximately 

L = G/Fo (2) 

as easily deduced from the tilted band picture. 

For the best understood case of a monochromatic driving F{t) — Fq — 
Fi cos{(Jit) where the driving frequency is in resonance with the Bloch fre- 
quency, ub = nuji, n = 1,2,..., a suppression of transport, denoted as a 
dynamic localization is found if the condition 

JnidFi/fkoi) = (3) 

is satisfied [1,4,5], where J„(-2) is the (ordinary) Bessel function. Otherwise 
transport with a velocity proportional to Jn{dFi/hu) is found. 

For the general case of a bi- or polychromatic periodic driving the particle 
dynamics is found to be quite complicated and unexpected effects have been 
observed. These effects may be used to manipulate the quantum dynamics. 
The aim of the present study is to investigate the transport and localization 
properties for a bi-chromatic driving extending previous studies by Liu and 
Zhu [6], Suqing et al. [7] and Yashima et al. [8]. In addition, we will briefly 
consider the even less explored realm of polychromatic driving 

We base our analysis on a Lie-algebraic approach introduced by one of 
the authors [9] (see also [10] for an extension to two space dimensions). 
Following Ref. [9], we introduce the abbreviations ft — dF{t)/h and g — 
—G/Ah which may also depend on time in general. Then the Hamiltonian 
can be conveniently expressed in terms of the hermitian position operator 
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N = T.e^\^)i^\ and the unitary shift operator K = + 1\: 

^H = g{K + K^)+ftN. (4) 
The dynamics of expectation values depends on the initial state 

mt = o)) = j2c,\i) (5) 

e 

characterized by the coherence parameters 

K = ^c|_iQ = lirle^'^ (6) 

L = ^c|_2Q = iLje^'' (7) 

e 

J = J2i2i-l)cUci=\J\e'r (8) 



In the following we will assume a symmetric normalized Gaussian state q ~ 
Q-f- /4o- ^hich allows for a broad initial distribution (a ^ 1) an approximate 
evaluation of the coherence parameters by replacing sums by integrals: 

K ^ e-^/«"' , L ^ e-i/2^' , J ^ . (9) 

The time evolution operator can be written as a product of exponentials [9] 

f/(t) = e"''"^ e"^^'^ e"^^*^^ , (10) 

with ^ ^ 

Vt= [ frdr and Xt = 9 [ e''^^ dr = \xt\ e-"^\ (11) 
Jo Jo 

An advantage of the product form ()10|) is the simple evaluation of matrix ele- 
ments and expectation values, as for instance the mean value of the position 

N, 

{N)t = {N)o + i{xt{K)o-X*t{K%) = {N)o + 2\K\\xt\ sin(0i - /c) (12) 

with the coherence parameter K = \K\e^'^ specified in For a symmetric 
initial state we have k = and an elementary calculation provides an upper 
bound of the mean transport velocity 

..an. = i^^^i^i^ < 2\gK\ . (13) 
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The time evolution of the width of the wave packet is given by [9] 

A%it) = A^(O) + 2|xt|'{l-|L|cos(20t-z/)-2|ir|2sin2(0,-/c)} 

+2\xt\{2{N)o\K\ sin(20t - k) + \J\ sm\(j)t - /i) } . (14) 

which simphfies to 

AUt) = A^(0) + 2\xt\^{l - \L\ cos(20O - 2\K\^ sin2(0,) } (15) 

for a real symmetric Gaussian initial wave packet. With \K\ = e''^/^"^ and 
\L\ = e~^/^'^ the dispersion depends on a as A^(t) — A^(0) ~ for 
large values of a. This leading order term can, however, be suppressed at 
times t for which the function Xt is purely imaginary, i.e. cos(20t) = —1 and 
sin^(0t) = 1. One can easily show that this implies 

AUt) - AUO) = \xt\' + 0{a-') (16) 

and the dispersion is strongly reduced at these times. This condition will be 
satisfied for the example of the shuttling transport discussed below in section 

o 



2 Time periodic driving 

In the following, we will assume a time-independent nearest-neighbor cou- 
pling g = g > and a periodic driving term 

ft = fo + ft , ft+T = ft , (17) 

where /o < is constant and ft vanishes on the average. This implies a 
Bloch frequency cub = /o- 

It can easily be shown by Fourier expansion that both functions in flllll 
can be separated into a linear growing and a periodically oscillating part 
(cf. [9]): 

Vt = fot + Vt, Xt = lt/2 + Xt (18) 

where 7 is different from zero if the driving period T is in resonance with the 
Bloch period Tb = 2tx/uj-b, 

T = nTB, neN. (19) 
The coefficient 7 controlling the transport properties is given by 

T 

^-inu^t-ivt^^ (20) 
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with uj = 27c/T. At time T the time evolution operator (fTUj) commutes with 
K and the simultaneous eigenstates of K and U{T) have eigenvalues e^'*'^ and 
e"'""^, respectively, where the quasienergy e depends on the quasimomentum 

e(K) = I7I cos^nd + arg(7) ) (21) 

denoted as the (quasi) dispersion relation. It may be of interest to realize that 
the possible values of 7 are bounded by the width of the dispersion relation 
for the non-biased system 

l7l<2|(7|, (22) 

a well known property of Fourier integrals. 

The mean value of the wave packet moves with an average velocity 

_ {N)Td 
•'trans rp ^trans 

-^trans = ^ Xt (^)o " Xt (^"^)o ) = - ItI slu (kc? + arg(7) ) , (23) 

where the initial condition is written as {K)o = e""^. As expected, this mean 
velocity can also be expressed as ftrans = de/dKl^ . Note that i^trans is bounded 
from above 

Wtrans < \l\ < '^^IqI OT Varans < 2d|5f| . (24) 

If the system parameters are adjusted such that 7 = 0, the band (PT|) 
collapses and the non-oscillatory term in xr vanishes with the consequence 
that there is no transport, i.e. we have dynamic localization. 



3 Bichromatic driving 

Let us consider a combined dc- and bichromatic ac-driving 

ft = fo~ fi cos(t^it) - /2 cos(u;2i + 5) (25) 

with a phase shift 6. Note that this function is aperiodic for incommensu- 
rable frequencies ui and uj2, which is not necessarily true for a a sum of two 
arbitrary periodic functions with incommensurable periods; an example can 
be found in [11]. For commensurable frequencies ft is periodic. 
With /o = ub and u = fi/ui, v = f2/uj2 'Vi'e obtain 

t 

rjt = J frdr = tOBt — u sin(co'it) — v sm(uj2't + 5) . (26) 




5 



The function Xt in eq. (fTT|) which determines the dynamics can be expressed 
in terms of Bessel functions using the generating function 



oo 

-,ix sm z 



E JA^)^'''- (27) 



e ' 

/i=— oo 

With the abbreviation cj^^^^ = uj-q ~ ^uji — 1^002 this yields 
t t 





f 

00 00 „ 

= E ^MW^-Me'^'y ^^e-"'""^. (28) 



The remaining integration is elementary, however one has to distinguish the 
cases of resonant and non-resonant driving. 
In the resonant case, we have 

ujb = ^uji + i>uj2 (29) 

for certain integers /x, z/. Let us denote the set of all indices /x, u satisfying 
(f^ by B. We then have 

Xt=7t/2 + 2g V J^(«)J,(^;)e-^^e-'--*/^in^^, (30) 

with 

7 = ^ = 2^7 JMUv)^'"'. (31) 

{/i,v)eB 

If B 7^ the first term in ()30p is growing linearly in time and thus dominating 
the long time dynamics. In this case we can observe transport. 

It is instructive to have a brief look at the monochromatic driving first. 
This case is recovered for t> = 0, observing Ji/(0) = 5^^ and hence 

Xt = lt/2 + 2gY^ J^(n)^ e--*/2 sin (32) 

with ijj^ = UB—fiuJi (cf. [1,9]) If the resonance condition ujb = fiuji is satisfied, 
we have 

7 = 2gJ^iu) (33) 
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and therefore an average transport with velocity proportional to the Bessel 
function J^{u) = J^^fi/ui). This linear transport is superimposed by an 
oscillating part, which is periodic with the Bloch period Tb for resonant 
driving. By adjusting the value of the driving amplitude /i one can change 
the direction of the transport or bring it to a standstill for J^(/i/co'i) = 0, 
the dynamic localization. 

For bichromatic driving, the situation is more complicated. We distin- 
guish two cases: 

Incommensurable driving frequencies ui and 002: 

If the ratio 002/ oji is irrational, the resonance condition ()29|1 leads to uo-b/ooi = 
fi + 1/002/(^1. This can only be satisfied for a special irrational value of the 
ratio ub/uji, for example by an appropriate dc- field Fq. These cases are of 
measure zero and typically there is no transport. If, however, the system pa- 
rameters are tuned to special values satisfying the resonance condition (j^ . 
then there exists only a single pair (/i, i/) e B as can be easily checked and 
we have transport with velocity 

-f = 2gJ^{u)Mv)e''''. (34) 

Commensurable driving frequencies Ui and UJ2: 

If UJ2/0J1 = q/p for integers q and p with no common divisor and if the 
resonance condition ()29|1 is satisfied for a particular pair (M, A^) one can find 
all solutions by {fi, v) = (M — qk, N + pk) with k E i.e. 

M = {{M -qk,N + pk),k eZ} . (35) 

This can be shown as follows: The resonance condition (j^^ can be rewritten 
as the Diophantine equation 

n = pji + qu with n = poj^/oJi E Ij . (36) 

A solution (/i, v) = (M, A^) of this equation always exists if p and q have 
no common divisor and can be found systematically by, e.g., the Euclid 
algorithm [12]. We therefore have an infinite number of solutions. Note that 
the choice of the special solution (M, A^) is arbitrary. For rational frequency 
ratios the function is determined by three integers, p, q, n, and we can 
express it in the convenient form 

7 = 2ge'^^ JP'%u,v] e'P^) (37) 

in terms of the two-dimensional, one-parameter Bessel function 

00 

JP^''{u,v;z) = J2 JM-qk{u)JN+pk{v)z'' (38) 
fc=— 00 
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ijj^jiox irrational 


ujT^juJx rational 




localization 


localization 


irrational 


(transport possible) 




localization 


transport 


rational 


(localization possible) 



Table 1: Transport properties for a bichromatically driven biased lattice 
with Bloch frequency uj-q. 



where {M,N) is an arbitrary solution of the Diophantine equation An 
introduction into the properties of two-dimensional Bessel functions can be 
found in [13]. 

The two-dimensional generalized Bessel functions fl38|) satisfy the inequal- 
ity \J^'''{u,v; z)\ < for n 7^ 0. This implies that the maximum trans- 
port velocity 

t^trans.max = ItI = 2|(7| | («, e'^') | (39) 

is bounded as I7 \<V2g<2g. 

For commensurable driving frequencies ui and U2, UJ2/UJ1 = q/p, and 
driving amplitudes /i /2 we have the following behavior: 

If the ratio ujb/^i is irrational, the condition ()29|) cannot be satisfied 
and we have no transport. 

If cub / 1^1 is rational with n = puj^ / Ui we typically have transport with 
transport velocity proportional to J^''^{fi/LJi, f2/uj2). This transport 
can, however, be stopped for special values of the system parameters 
and dynamic localization is observed. It should be noted that this case 
is always met for pure ac-driving (cub = 0). 

The transport and localization properties of a bichromatically driven system 
are summarized in Tabled The following section considers some exemplary 
cases in more detail. 

4 Case studies for resonant driving 

In this section, we will discuss the transport and localization properties in 
some more detail for the most interesting case of resonant driving ub = fiu!i + 
VUI2 ■ The driving frequencies are assumed to be commensurable, UJ2/UJ1 = p/q 
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Figure 1: Transport coefficient 7 as a function of the parameters u = fi/uji 
and V = f2/^2 for a bichromatic driving with UJ2 = 2uJi at resonance with 
the Bloch frequency u^, i.e. ui + UJ2 = uub- Shown are examples for n = 1 
(left and middle) and n = 29 (right), where the figure in the middle shows 7 
for a fixed parameter v = 1. 

with coprime integers p and q, where, w.l.o.g., we can choose p as the smaller 
one of these numbers. Let us recall that the ac-driving period T is equal 
to an integer number of Bloch periods, T = nT^, where n can be expressed 
as n = pM + qN with integer (M, A^) G B. We will also assume that the 
force term is symmetric in time, /(— t) = f{t), i.e. the phase shift 6 in ()25|1 is 
equal to zero, with the consequence that the two-dimensional one-parameter 
Bessel functions reduce to simple two-dimensional ones [13]. The transport 
properties of the system are then determined by the parameter 

7 = 2gjr{u,v). (40) 

with variables u = fi/ui and f = /2/CU2. In this case 7 is real and its sign, 
the sign of the Bessel function, directly gives the direction of transport. Here 
we are in particular interested in the parameter values leading to dynamic 
localization. These values are given by the zeros of J^''^{u,v), i.e. the nodal 
lines of the two-dimensional Bessel function. 

In previous studies the special case of the bichromatic driving with 002 = 
2uJi attracted most attention and we will start with this special case in the 
following. For simplicity we fix = 1 and ui = 1. The resonance condition 
for the static field (jHUj) then reads n = ub/uji = fo/ui. In our first examples 
we fix the amplitude of one driving as /2 = UJ2, i-e v = 1, and consider the 
case lvb = uji, i.e. n = 1. The transport coefficient given by 7 = 2j|'^(-u,t>) 
is shown in figure Q in dependence of u and v (left-hand side) and for v = 1 
fixed (middle). Numerically we find a global maximum at u = —4.68 and a 
nodal line at u = —6.49 for v = 1. 
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Figure 2: Left: Dynamics of a gaussian wave packet with initial momentum 
K,Q = — 7r/2 for a bichromatic driving with amphtudes v = 1 and u = —4.68 
in a greyscale plot. Right: Position expectation value {N(t)) for kq = —tt/2 
(solid line) and = (dash-dotted line) in comparison with the transport 
coefficient jt (dashed line). 

We consider the time evolution of a broad Gaussian wavepacket 

m = o))r^Yl H'/i^^y + M 1^) (41) 

e 

with a = 10. For such a broad gaussian wavepacket the coherence parameter 
K is given by ~ 1 and k, ^ kq- The left-hand side of figure |21 shows the 
dynamics of such a wave packet for kq = — 7r/2, which is rapidly transported. 
The right-hand side shows the expectation values {N{t)) of the position for 
Kq = and Kq = — 7r/2. The transport coefficient 7t is shown as a dashed 
line for comparison. The sine in equation (|T2|l is approximately one for 

= — vr/2 such that the fastest possible transport is observed which is given 
by jt. For = the sine in equation (fT^ is approximately zero so that 
the wavepacket stays at rest. However, this localization is only due to the 
special initial state chosen in this example. Depending on the value of Kq, 
the wavepacket will be transported with a velocity in the interval [—7,7]. 
Localization, i.e. ftrans = is only found for kq = 0. 

The situation is different for u = —6.49, for which jI''^{u,v) = 0, i.e. 
dynamical localization is found. Now every wavepacket will be localized 
regardless of its initial momentum k. This is illustrated in figure El on the 
left where the position expectation value {N(t)) is shown for kq = and 
Ko = —7t/2. The right hand side of figure El shows the time evolution of 
the width of the wavepacket with kq = for u = —4.68 (transport) and 
u = —6.49 (dynamical localization). As the dispersion is also governed by the 
function \xt\, no broadening is observed in the case of dynamical localization. 
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Figure 3: Left: Dynamical Localization. For u = —6.49 any wave packet is 
localized around n = regardless of the initial momentum. The evolution 
of the position expectation {N{t)) is shown for kq = —tt/2 (solid line) and 
Ko = (dash-dotted line). Right: Dispersion of a gaussian wave packet with 
Ko = —7t/2 for u = —4.68 (transport, solid line) and u = —6.49 (dynamical 
localization, dashed). 

Next we consider an example, where the field amplitudes are large in 
comparison to the driving frequencies. Then n = /o/u; ^ 1 as well as 
u, f -C 1 and asymptotic approximations for the relevant Bessel functions 
are available [13]. In fact we consider a driving amplitude v = —20 and a 
static field strength with n = 29. Dynamical localization is observed if the 
driving amplitudes u is chosen such that jI''^9{u,v) = 0. Explicit estimates 
for these values are derived from the asymptotic approximations of this Bessel 
function. For n < 2\v\ one finds 

u.l---= I + l)f n even . = g ±1 ±2 (42) 

V 2 Av \ JTT n odd ' ^ ^' ■ ■ ■ ^^^^ 

and for n > 2\v\ with v < 

\ 7T 77" 

"V2"4^= (""^^-^'^ ' ^■ = 0'±1'±2,.... (43) 

The smallest non-zero values of u for which dynamical localization is pre- 
dicted by these formula predict are ui = 3.38 and U2 = 6.77, while we find 
Ui = 3.37 and U2 = 6.75 numerically. Figure 0] shows the time evolution of 
the position expectation value for ui = 3.37 (left) and for M3 = 5 (right), 
each for Kq = — 7r/2. Clearly the wave packet stays and rest and does not 
show any systematic dispersion for ui = 3.37. Transport is found for ^3 = 5, 
where the mean transport is given by '-ft. However, transport is much slower 
than for n = 1 as illustrated in figure El 
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Figure 4: Dynamical Localization for n = 29, v = —20 and u = 3.37 (left) 
and transport for -u = 5 (right). The mean transport is given by 'jt, which is 
plotted as dashed line for comparison. 

5 Polychromatic driving 

In view of the rich transport and localization behavior for bichromatic driv- 
ing, one may expect additional surprising effects for general time periodic, 
however polychromatic driving. This is indeed the case, however, only few 
such studies have been reported up to now in particular for a periodic rect- 
angular force [2,14]. This case allowing a closed form solution, is discussed 
in the following section. The general case is briefly outlined in sect. 15.21 

5.1 Flipped fields and shuttling transport 

Let us consider a rectangular force term 

= I /, aT<t<T C"" 

with < a < 1 and f{t + T) = f{t). For this flipped static force excitation 
the time-averaged field is 

fo = ah + hh , 6 = l-a, (45) 

where we can assume /i > 0. For the choice /2 = — /i, the force (jl^ reduces 
to the case studied in [14] (note that there the time scale is shifted). The 
even more specialized case a = 1/2 with period T = 27r//i, the Bloch period 
for a constant field /i, is discussed in [15] and briefly at the end of [2] (note 
that in this case we have /o = 0). 
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For the flipped field excitation (jUJ, the functions rjt and Xt in (fTTj) can 
be calculated by elementary integration: 



fit 0<t<aT 
^* " ^ fiaT + /2(t -aT) aT <t<T ^^^^ 

rjt+T = Vt + VT , Vt = fioT + /26T = /ot . (47) 



t|- (1 - e-'^^*) < t < aT 

= ' ' -i/,aT (48) 

A f 1 _ e-'^^"^) + ^-^-^ ( 1 - e-^^^(*-"^)) aT<t<T 
1/1 V ^ 1/2 V / 

and 

Xt+T = e-'^^^Xt + XT (49) 

XT = Tf (1 - e-^-^) + :f (e-/-^ - e^/"^) . (50) 
1/1 ^ ^ 1/2 V / 

With /o = cub and u = 27r/T, we see that for the case of resonant driving, 
ub = nu, we have /qT = 27m , XkT = kxT for k = 1, 2, . . .. This implies an 
overall transport of the expectation value {N)t with a velocity determined 
by 

For a/iT = 2/i7r with /i = 1, 2, ... we have 7 = 0, i.e. dynamical local- 
ization. The maximum transport velocity 



^trans | T I 



4(7/1 1\ . a/iT 

T (/:-/•) ^^^^ 



< 2|(?| (52) 



is obtained when the phase is adjusted as nd + arg(7) = nd — afiT/2 = 
(2/i + l)7r/2 according to (1^ . 

These findings can be easily explained in terms of Bloch oscillations for 
the time intervals with constant force /i and /2 with Bloch periods Tbi = 
27r//i and Tb2 = 27r/|/2|, respectively. If a/iT = 2/i7r with /i G N, we have 
^1/2!^ = 2i^7r with z/ G N for resonant driving, /q = af\T + bf2T = 2mT. The 
wave packet carries out fi full Bloch oscillations in the first time interval and 
u in the second and there is no no directed transport. For a/iT = (2/i + l)7r, 
however, we also have a\f2\T = (2z/ + l)7r [fi, z/ G N) and therefore a motion 
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Figure 5: Numerical time propagation of a Gaussian wave packet in a cos- 
potential for a flipped field fi = —f2 = 0.0003 and a = 0.5. 



over a distance Li = A/Fi = Agd/fi in the first time interval aT and L2 = 
A/F2 = 4gd/f2 in second time interval bT, however in the opposite direction. 
The average transport velocity is 



_ L1 + L2 _ 4gd 

•'trans rp rp 



1 1 
Ji~T2 



(53) 



in agreement with 

Let us briefly consider the most simple case where a static field / with 
Bloch period Tb = 27r// is flipped after each half of a Bloch period, i.e. we 
have /i = — /2 = /, a = 1/2, /o = and a period T = Tb. This yields a 
transport velocity 

^^trans = — , (54) 
71 

which is by a factor of 2/7r smaller than the maximum possible value 2\g\. It 
is remarkable that the velocity of this field induced transport is independent 
of the field amplitude /. 

The analysis above is based on a tight-binding model. Numerically exact 
wave packet propagation can be used to check the validity of these predictions 
for realistic potentials. As an illustrating example, figs. El and IHl show the 
propagation of an initially broad Gaussian wave packet in position space, 

^{x,t = 0) ~ exp(-a;V(2s)^) , s = 157r (55) 
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Figure 6: Numerical time propagation of a Gaussian wave packet for a 
flipped field /i = -/s = 0.0003 and a = 1/4. 

(we use units with ^ = 1) in a potential V{x) = (1/8) cosa; for a flipped 
F{t) field with \F\ = 0.0003 and /i = —f2- In the first case shown in fig. El 
the static force is flipped after each half-period of the Bloch time (a = 1/2), 
i.e. the average field /o vanishes (see also [2] for a similar calculation). One 
observes that the wave packet moves almost dispersionless with a velocity 
Krans = "78 (i/Te, which agrees with the prediction of the tight-binding 
model for the actual bandwidth A = 0.0741. Note that the direction of 
transport is determined by the direction of the Bloch oscillation at time t = 0, 
i.e. in the negative direction in the present case with initially positive /. 
In order to demonstrate this clearly, fig. IHl shows the same system as before, 
however for a doubled excitation period T and a = 1/4, i.e. a negative average 
force /q. The wave packet moves in the same direction as before, i.e. opposite 
to the average gradient of the potential. Note that the transport velocity is 
reduced by a factor of 2, because of the increased value of T. This behavior 
can be easily understood by observing that in the second half of the excitation 
period, the system undergoes a full Bloch oscillation without a net transport. 

5.2 General polychromatic driving 

The two simple cases analyzed above, mono- or bichromatic fields or binary- 
flipped piecewise constant fields allow a quite extensive analytical analysis. 
Let us finally briefly address a general polychromatic driving field following 
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[9]. For simplicity we will assume f{t) = f{t + T) to be symmetric in time, 
/(— t) = f{t), with Fourier expansion 

oo 

ft = fo-^fm cos{mujt) , UJ = 27r/T . (56) 

m=l 

Let us recall that in the units chosen here the dc-component is equal to the 
Bloch frequency /q = cjb, and the function rjt in eq. (fTTj) is 



Vt 

This yields 

Xt = 



UBt -^Pm sm{mut) , Prn = fm/muj . (57) 



m=l 



Jo Jo 

+00 „t 



{51 



in terms of the infinite- variable Bessel functions [16, 17] 

oo +00 

exp (i ^ Pm sin mn) = ^ J.({/3,r^}) e"^'^" . (59) 

m=l i/=— 00 

For resonant driving, uj-q = nu, the final result is 

Xt = gU{Pn.})t + 2gy2 Jumm}) - e-'-^*/^ sin ^ , (60) 

where we have set = ub — voj. In the non-resonant case the first term is 
absent and the sum extends over all integers. Specializing again to mono- or 
bichromatic driving this result reduces to the ones derived in sect. 01 

For resonant driving, ujb = nuj, we have uj,^ = [n — v)uj and the oscillating 
part in ()60|) vanishes at times which are integer multiples of the driving period 
T. The average transport velocity is then given by 

t^trans = ItI = = \g Jn{{(3m})\ ■ (61) 

Dynamical localization effects will be observed if the system parameters are 
tuned to a zero of the infinite-order Bessel function Jn{,{Pm\)- 

The closed form result for the transport velocity in ()61|) is a generalization 
of the formula ()39|) derived for bichromatic driving. For the very special case 
of a flipped rectangular field, the closed form solution presented in section 
15.11 is. of course, superior to the general solution. 
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6 Concluding remarks 



We have shown, that the dynamics of a quantum particle in a one- dimensional 
periodic potential (lattice) under the action of a static and time-periodic 
field shows quite intricate transport and localization phenomena, even in the 
relatively simple case of bichromatic driving which depends sensitively on 
the communicability relations between the two excitation frequencies and the 
Bloch frequency. The theoretical analysis based on a nearest-neighbor tight- 
binding model allows a convenient closed form description of the transport 
properties in terms of generalized Bessel functions. In particular, conditions 
for the system parameters leading to localization are derived. 

The case of polychromatic driving is also discussed, in particular for 
flipped static fields, i.e. rectangular pulses, which can support an almost 
dispersionless transport with a velocity independent of the field amplitude. 

Such fields, which can be quite easily realized experimentally, offer inter- 
esting possibilities for the control of quantum transport and the manipulation 
of wave packets as for instance cold atoms in optical lattices. 

Finally it should be noted that it offers also some advantage to modify 
the spatial period of the lattice. For example, a double periodic lattice, 
a superposition of two periodic structures with period d and 2d has been 
studied recently [18,19]. It allows controlled Landau- Zener tunneling which 
can also be employed for a manipulation of matter waves. 
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